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to ' ABSTRACT 

c3 : 

■ Context. Lines of diatomic molecules are typically much more temperature and pressure sensitive than atomic lines, which makes them ideal, 
I complementary tools for studying cool stellar atmospheres as well as the internal structure of sunspots and starspots. The FeH F 4 A-X 4 A system 

■ represents such an example that exhibits in addition a large magnetic field sensitivity. However, the current theoretical descriptions of these 
QO ' transitions including the molecular constants involved are only based on intensity measurements because polarimetric observations have not 

i been available so far, which limits their diagnostic value. Furthermore, the theory was optimized to reproduce energy levels and line strengths 
£NJ ■ without taking the magnetic sensitivities into account. 

CiD Aims. We present for the first time spectropolarimetric observations of the FeH F 4 A-X 4 A system measured in sunspots to investigate their 
. diagnostic capabilities for probing solar and stellar magnetic fields. In particular, we investigate whether the current theoretical model of FeH 

■ can reproduce the observed Stokes profiles including their magnetic properties. 

[ — ' Methods. The polarimetric observations of the FeH F 4 A-X 4 A system in Stokes / and V are compared with synthetic Stokes profiles modeled 
C^> with radiative transfer calculations. This allows us to infer the temperature and the magnetic field strength of the observed sunspots. 

' Results. We find that the current theory successfully reproduces the magnetic properties of a large number of lines in the FeH F 4 A-X 4 A system. 
, — In a few cases the observations indicate a larger Zeeman splitting than predicted by the theory. There, our observations have provided additional 

■ constraints, which allowed us to determine empirical molecular constants. 

<-h ' Conclusions. The FeH F 4 A-X 4 A system is found to be a very sensitive magnetic diagnostic tool. Polarimetric data of these lines, in contrast 
to intensity measurements, provide us with more direct and detailed information to study the coolest parts of sunspot and starspot umbrae, and 
cool active dwarfs. 
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1. Introduction 

Many lines of diatomic molecules are not only good tempera- 
ture and pressure indicators but also excellent magnet i c field 
diag nostics, as has been s h own b y iBerdvugina et al. ( 2000b 
and IBerdvugina & Solanki (2002) in their overview of mag- 
netic properties of molecules observed in spectra of sunspots 
and cool stars. 

Iron hydride (FeH) is one of the most sensitive indicators 
of magnetic fields in cool atmospheres, but s o far observa- 
tions h ave been restricted to intensity spectra. IWallace et al. 
dl998l) revealed the remarkable magnetic sensitivity of many 
lines of the FeH F 4 A-X 4 A system in sunspot observations. 
Valenti et al.l (12001 ) detected Zeeman broadening of FeH lines 
in an active M dwarf and illustrated the usefulness of these 



lines in studying stellar magnetic fields. Because of the lack 
of a theoretical description of the FeH molecule, they modeled 
the stellar spectrum employing the sunspot measurements by 



Wall ace et al.l (119981) . First modeling of synthetic Stokes pro- 



files of FeH lines indicated that this particular electronic sys- 
tem will be a powerful tool for diagnosing solar and stellar 
magnetism once the spin -coupling constants become available 
dBerdyugina et al.ll2003h . 



The set of mis sing molecular constants was provided by 



Duli ck et al.l P003). They determined the spin-orbit (A) and 
spin-spin (A) constants for th e FeH F 4 A-X 4 A sys tem, using the 



laboratory measurements by Phil lips et al.l (119871) . They found 
that the spin splittings for the observed vibrational levels could 
be reproduced within an error of 0.5 cirT 1 by a single pair of A 
and A values for each electronic state. However, the analysis by 
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Dulick et al.l Q2003) aimed at modeling energy levels and line 
strengths without considering magnetic characteristics. Only if 
the theory correctly describes the magnetic properties we can 
exploit the full diagnostic value of FeH. 

Here we present first polarimetric measurements of the FeH 
F 4 A-X 4 A system which provide novel constraints and tests for 
the theory. They were recorded in a sunspot and are described 
in Sect. [2] We carried out a perturbation analysis for the FeH 
F 4 A-X 4 A system to account for the Zeeman effect as described 
in Sect. [3] Molecular Stokes line profiles in the presence of 
magnetic fields are calculated and compared with the sunspot 
observations in Sect. |4] Our results are based on the molecu- 
lar constants determined by lDulick et al.l d2003l) to test the cur- 
rent theory. We conclude that polarization measurements of the 
considered FeH system and their successful theoretical inter- 
pretation strongly constrain molecular parameters and greatly 
advance our ability to explore magnetic fields in cool stellar 
atmospheres. 

2. Observations 

Observations were performed at Istituto Ricerche Solari 
Locarno (IRSOL) with the 45 cm aperture Gregory Coude 
telescope, the high spectral resolution Czerny Turner spectro- 
graph, and a two beam exchange polarimeter. On 23rd August 
2004 registrations were done in a sunspot (NOAA 0663). 
Wavelength intervals of about 5 A were observed around 9902 
A, 9906 A, 9976 A, 9979 A, 10059 A, and 10062 A. 



The two beam exchange technique by ISemel et al. I (1 19931) 
had already been successfully applied at IRSOL f or polarimet- 



ric ob servations in the blue part of the spectrum (Bianda et al 



1998). To observe the FeH lines in the near infrared we used the 
origin al Semel polarimeter described by lSemel & Lopez Ariste 
(1200 ll) without collimation lenses. 

The polarization analyzer consists of a set of motorized 
Fresnel Rhomb retarders installed in front of a Savart Plate 
beam splitter. The retarders are rotated in order to separate the 
incoming beam into two beams corresponding to I + X and 
I - X with X being Stokes Q, U, or V. Thus, two images of 
opposite polarization of the same solar region were projected 
on the spectrograph slit jaws. To obtain a single Stokes mea- 
surement two such spectral recordings were registered with a 
CCD subsequently changing between the polarization state of 
the two beams with a suitable rotation of the retarders. On the 
spectrograph slit, two images with orthogonal polarization are 
projected, each subtending 35" in the spatial direction along 
the slit direction. The position on the spectrograph slit was cho- 
sen in order to maximize the length of the segment covering the 
sunspot umbra. The profiles in Fig. 1-3 are obtained averaging 
the coolest part of the umbra spectrum over 7". The exposure 
time for each frame was 5 seconds. A complete cycle of Stokes 
vector observation required six frames and took about 2 min- 
utes, leading to a spatial resolution of 2 arcsec. Two cycles were 
performed for each observed spectral interval. 

Dark frames were recorded at the beginning and at the end 
of the observation run. Flat field observations were performed 
near the solar disk center on quiet re gions. Details of the data 
reduction technique are described by iBianda et alj (1998). To 



correct for the instrumental polarization the polarization val- 
ues obtained in the flat field registrations at disk center were 
subtracted from those measured in the sunspots. 

3. Theoretical approach 

If a molecule possesses a nonzero magnetic moment, it inter- 
acts with an external magnetic field and causes a precession 
of the total angular momentum J about the field direction. As 
a result one can observe the Zeeman effect, i.e. the magnetic 
splitting of energy levels in the molecule. 

The energies of magnetic components depend on how the 
angular momenta are coupled to the rotation of the molecule. 
The FeH F 4 A-X 4 A system is produced by transitions between 
two electronic states with couplings of the angular momenta 
that are intermediate between the limiting Hund's cases (a) and 
(b). The perturbation calculation of the molecular Zeeman ef- 
fect is based on the effective Hamiltonian for this transition, 
which can be represented by two parts, < H ab and "H H , 



K = <H ah + H H 



(1) 



the former performing the transformation of the Hund's case 
(a) wavefunctions into the intermediate case wavefunctions and 
the latter describing the interaction with the external magnetic 
field. The matrix elements of the Hamiltonian < H ab include the 
spin-orbital, spin-spin, rotational, and centrifugal distortion in- 
teractions in the electronic st ate described with H und's case (a) 
wavefunctions as a basis set ( Duli ck et al ] |2003l) . By diagonal- 
ization of the total Hamiltonian "7Y we obtain the energy levels, 
the transition strengths, and t he Zeeman pattern of th e magnet- 
ically perturbed system (see lBerdyugina et al ] l2005l) . 



4. Stokes profiles in the Zeeman regime 

By implementing the Hamiltonian given in Eq. £0, we can now 
calculate molecular line profiles in the presence of magnetic 
fields in th e Zeeman regime, using the radiativ e transfer code 
STOPRO dSolankil Il987t ISolanki et al.l Il992l iFrutiger et al ' 



2000; Berdvu ginaetal.ll2003l l2005h . These synthetic Stokes 



profiles are compared with profiles observed in sunspots. For 
this purpose we vary only the model atmosphere and the mag- 
netic field in our fitting procedure . Atmospheric models are 
taken from lHauschildt et al.l ( 1999b . To account for the pho- 
tospheric stray light, we combined calculated sunspot and pho- 

tospheric spectra using various spot fillin g factors. 

With the molecular constants from Dulick et al. (2003) 



many spectral lines with various A/, J, and D. can be repro- 
duced, J being the total angular momentum and Q the quan- 
tum number for the component of the total electronic angular 
momentum along the internuclear axis of the molecule. 

In Figs. 1 and 2 we present such examples of observed FeH 
Stokes / and V profiles compared wit h the calcula t ed spe ctra. 



(1998) and 



Stokes / observations are taken from lWallace et al. 
IRSOL (upper panels), whereas Stokes V spectra (lower pan- 
els) were obtained only and for the first time in observations at 
IRSOL. The fact that the spots observed at IRSOL were neither 
very large nor cold is reflected in the shallow depth of the FeH 
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Fig. 1. FeH Stokes / and Stokes V line profiles (upper and lower 
panel, respectively). Solid lines represent synthetic Stokes pro- 
files; observed Stokes / and V profiles are shown with dotted 
lines . The top Stokes / s pectrum shows a sunspot observation 
from lWallace et aL ( 19981) . while the bottom Stokes / spectrum 
(shifted for the sake of clarity) as well as the Stokes V spec- 
trum correspond to our own sunspot observatio ns at IRSOL. 



Synthetic profiles are based on the theory by Dulick et al 



(2003) and were calculated with a field strength of 3 kG, um- 
bral mod els with effective te mperatures T e ff of 4300 K and 
4800 K dHauschildt et alJll999h . and filling factors of 0.95 and 
0.8 to fit the Wallace and IRSOL data, respectively. 



intensity profiles as compare d to the spectra from the sunspot 
atlas by I Wallace et al.l 11998). This leads to smaller filling fac- 
tors in the IRSOL spectra since the influence of stray light 
from the photosphere becomes more relevant. Nonetheless, the 
Stokes V profiles of the FeH lines have very clear, easily mea- 
surable signals with a degree of circular polarization of the or- 
der of 1 % to 3%. The magnetic field strengt h was foun d to have 
the same value of 3 kG in both the Wal lace et al.l (119981) and the 
IRSOL data sets. In both cases the magnetic field vector was 
assumed to be aligned along the line of sight in our modeling, 
pointing away from the observer, although this direction can of 
course only be determined for the IRSOL data thanks to the 
spectropolarimetric measurements. Since FeH lines effectively 
form at log T(lfi) = -2.5, according to our model, the inferred 
magnetic field corresponds to this depth. 

The main features of both Stokes / and V could be re- 
produced, suggesting that our calculations of the molecular 
Zeeman spli t ting in these FeH lines based on the theory by 
Dulick et alJ (|2003) are realistic. Lines of the £>-branch with 
J < 4.5 exhibit particularly strong mag netic sensitivity as was 
expec ted from earlier investigations dBerdvugina & Sol anki 
20021) . Also the P and R lines with Q < 2.5 and J < 10.5, such 
as those presented in Figs. Q] and El can b e properly modeled 
with the molecular constants from 



Fig. 2. The same as Fig. Q]but for a different wavelength re- 
gion. Synthetic pro files (solid) are based on the theory by 
Dulick et al.1 ( 12003!) and were calculated with a field strength 
of 3 kG, umbral mod els with 7/ eff of 4000 K and 4500 K 
dHauschildt et al.lll999t) . and filling factors of 0.97 and 0.8 to 
fit the Wallace and IRSOL data, respectively. 



opposite sign of the effective Lande factors. This represents an 
example of important additional constraints for the theory that 
remain inaccessible by intensity measurements and can only be 
provided by spectropolarimetric observations. 

Some of the observed line profiles arising from the FeH 
F 4 A-X 4 A system cannot be modeled correctly when employ- 
ing the theory by iDulick et al.1 ( 2003 ). e.g. the R4 line£] in the 
region from 9900 A to 10000 A. Figure [3] illustrates examples 
of such R4 lines which would be ideally suited for magnetic 
field diagnostics thanks to their strong splitting and peculiar 
profile shapes. However, their effective Lande factors and thus 
their Zeeman splitting as calculated with the current theory are 
about one order of magnitude too small. These lines with a 
large LI = 3.5 and J > 10.5 seem to be affected by unknown 
perturbations neglected in the theory such that neither the inten- 
sity nor the polarization profiles are properly reproduced. For 
those transitions we searched for ideal spin-orbit interaction 
constants A with ax 2 minimization. Best fits to the Stokes pro- 
files were obtained for A" = -30cm _1 and A' = -330cm _I , 
for the lower and upper level respectively. 

Unknown perturbations appear to have a small influence on 
lines of the Q-branch (AJ = 0) in which the /-number of the 
lower and upper level coincide so that perturbations seemingly 
cancel out. For the P- and R -branches such cancelations do not 
occur and perturbations might become a major issue for lines 
with Q > 2.5 and J > 10.5, which will be investigated in a 
forthcoming paper. Such a behavior is similar to that d iscussed 
for the TiO y' system by lBerdyugina & Solan ki (2002). 



t lDulick et alJ (120031) . Note in 
particular the opposite sign of the Stokes V profiles in the P- 
and /Mines as compared with the g-li nes i which is due to the 



1 Here we denote rotational branches with indices 1, 2, 3, and 4 as 
corresponding to transitions between levels both having L~l = 1/2, 3/2, 
5/2, and 7/2. 
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Fig. 3. The same as Fig. [TJ for a different spectral window 
containing FeH lines that suffer from unknown perturbations. 
Observations (dotted) are compared with two different models, 
the first being based on our empirically determined spin-orbit 
constants (solid), and the se cond one (dashed) assuming the 
theory bv lDulick et al. (2003). Both synthetic spectra were cal- 
culated with a field s trength of 3 kG, umbral models with T e ff of 
4300 K and 4800 K dHauschildt et all 19991) . and filling factors 
of 0.98 and 0.9 to fit the Wallace and IRSOL data, respectively. 



5. Conclusions 

The first polarimetric observations of the FeH F 4 A-X 4 A system 
and their successful modeling open an opportunity to fully ex- 
ploit the diagnostic capabilities of this molecular system, which 
represents a powerful tool for solar and stellar magnetic field 
studies. 

We have described the synthesis of molecular Stokes pa- 
rameters and showed that it was possible to fit a wide range 
of available spectral lines in different wavelength regions with 
various A/, J, and £1 numbers with the set of molecular con- 
stants provided by IDulick et al.l (120031) . especially lines with 
small /-numbers. Our analysis has demonstrated that the the- 
ory produces not only the true Zeeman splitting but also the 
right sign for the Stokes V profiles. Since the theory for these 
lines seems to be understood, they represent good candidates 
for determining the thermal and magnetic structure of cool at- 
mospheres. 

Nevertheless, interesting spectral regions exist, such as the 
region around 9900 A , where a strong splitting is observed that 
cannot be modeled with the available constants. It seems that 
an unknown perturbation, which is not taken into account by 
the applied Hamiltonian, affects some of the magnetically very 
sensitive lines. However, using spin-orbit constants which were 
determined empirically, we were able to reproduce remarkably 
strongly split Stokes / and Stokes V profiles observed in a 
sunspot spectrum, such as the /?4(13.5), /?4(14.5), and /?4(15.5) 
lines. 



We conclude that our ability to model the FeH F 4 A-X 4 A 
system based on a set of theoretical and empirical molecular 
constants provides a new valuable source of information about 
the physical parameters in observed objects. 
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